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ABSTRACT 

The impact of the combination of rotation and magnetic fields on oscillations of stellar 
fluids is still not well known theoretically. It mixes Alfven and inertial waves. Neutron 
stars are a place where both effects may be at work. We wish to decipher the solution 
of this problem in the context of r-modes instability in neutron stars, as it appears 
when these modes are coupled to gravitational radiation. 

We consider a rotating spherical shell filled with a viscous fluid but of infinite 
electrical conductivity and analyze propagation of modal perturbations when a dipolar 
magnetic field is bathing the fluid layer. We perform an extensive numerical analysis 
and find that the m — 2 r-mode oscillation is influenced by the magnetic field when the 
Lehnert number (ratio of Alfven speed to rotation speed) exceeds a value proportional 
to the one- fourth power of the Ekman number (non-dimensional measure of viscosity) . 
This scaling is interpreted as the coincidence of the width of internal shear layers of 
inertial modes and the wavelength of the Alfven waves. Applied to the case of rotating 
magnetic neutron stars, we find that dipolar magnetic fields above 10 G are necessary 
to perturb the r-modes instability. 
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1 INTRODUCTION 

Numerous astrophysical systems exhibit a pulsating behav- 
ior that can be significantly influenced by both magnetic 
field and rotation. The rapidly oscillating Ap (roAp) stars, 
the magnetic white dwarf stars and neutron stars as well as 
planetary cores fall into the above category. 

In neutron stars, however, possible astrophysical impli- 
cations of r-modes instability have motivated extensive in- 
vestigations of this mode during the past few years (see An- 
dersson 2003 for a review), r-modes belong to the class of 
inertial modes that arise in rotating fluids due to the Cori- 
olis force. It was shown by Andersson (1998) and Friedman 
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& Morsink (1998) that these modes easily couple to gravita- 
tional radiation and become unstable, allowing the neutron 
stars to lose their angular momentum. 

The foregoing instability may however be weakened or 
even suppressed by all the dissipative mechanisms which 
couple to an r-mode oscillations. Much work has thus been 
devoted to the analysis of the various mechanisms which 
may damp r-modes. Recently, vortex-mediated mutual fric- 
tion of superfluids was investigated by Haskell et al. (2009) 
as well as the effects of hyperon bulk viscosity (Haskell & 
Andersson 2010), but both effects do not seem to be able 
to influence the instability at low temperatures. For these 
temperatures, the Ekman layer which forms below the crust 
of a neutron star, remains the most important source of 
dissipation (Bildsten & Ushomirsky 2000; Rieutord 2001; 
Glampedakis & Andersson 2006). 

However, it is well-known that the magnetic field is 
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an important component of a neutron star. It is therefore 
clear that fluid flows may be seriously influenced by this 
field and that other channels of dissipation for the r-mode 
instability may exist through this component. Much work 
was devoted to investigate the modifications induced by 
a magnetic field on these kind of modes. Rezzolla et al. 
(2000, 2001a,b), for instance, showed that strong magnetic 
field, beyond 10 10 G may weaken the r-mode instability suf- 
ficiently so as to make the generated gravitational waves 
undetectable. Mendell (2001) and Kinney & Mendell (2003) 
focused their study on the influence of the magnetic field 
on the Ekman layer flow. They conclude that magnetic field 
larger than 10 12 G completely suppress the instability. How- 
ever, Lee (2005), using a dipolar magnetic field covering the 
surface of a neutron star modeled as a N = 1-polytrope, con- 
cludes that much stronger magnetic fields, over 10 14 G, are 
necessary to suppress the instability through magnetic per- 
turbations. Let us note that beside the suppression of the 
foregoing instability, magnetic fields can also directly spin 
down a neutron star through the classical process of mag- 
netic braking (i.e. Ho & Lai 2000). More recently, Cuofano 
& Drago (2010) and Bonanno et al. (2011) investigated the 
role of magnetic field generation by the unstable r-mode, 
showing that a very strong toroidal magnetic field can be 
generated, which in turn can modify the instability. 

The aim of the present paper is to explore furthermore 
the channels of dissipation for the unstable m = 2 r-mode of 
a rotating neutron star, especially when viscosity and mag- 
netic fields are both present. Indeed, a possibility that has 
not been considered by previous studies (like the one of Lee 
2005), is the fact that an unstable r-mode of a spherical 
layer, might develop internal shear layers thanks to the mag- 
netic field action. Since fluid layers of different natures are 
expected due to phase transitions of nuclear matter, the ex- 
istence of internal shear layers are also very likely. We shall 
see that this leads to a threshold value around 10 14 G for 
magnetic fields to noticeably perturb the r-mode instability. 

Rotating fluid layers bathed by a magnetic field are 
however not specific to neutron stars. They are also found 
in planetary core, like in the Earth's core. Thus, in order to 
make the following study of more general interest, we shall 
consider a very simplified model of neutron star, neglecting 
relativistic or superconducting effects. We thus extend the 
results of Rieutord (2001) to the case where a dipolar mag- 
netic field perturb the fluid flow (in the limit of very large 
magnetic Prandtl numbers). This is meant to be a simple 
configuration where combined effects of magnetic fields and 
viscosity can be studied. 

The paper is organized as follows: In Sect. 2, we recall 
the basic physical ingredients of the model and then (sect. 3) 
briefly explain the numerical strategy. Numerical results for 
the m — 2 r-mode coupled with Alfven waves are discussed 
in Sect. 4. Conclusions follow. 



2 THE MODEL 

We consider a rotating star modeled as an infinitely electri- 
cally conducting core surrounded by a spherical layer of fluid 
itself limited by an outer solid crust. The ratio of the inner 
core radius to the outer radius R is rj. The kinematic vis- 
cosity and magnetic diffusivity of the fluid are respectively 
v and v m . 

The star is rotating with uniform angular frequency 
ft = Qe z along the z-axis. The core is supposed to be the 
source of a permanent axisymmetric magnetic dipole cov- 
ering the whole layer. The symmetry axis of the magnetic 
field is along the z-axis. Note that such a partition of the 
star is necessary in order to avoid the problem of the defini- 
tion of the magnetic field in the core of the star, the dipole 
field being singular at the center. The fluid in the layer is 
taken to be incompressible; we thus eliminate phenomena 
related to compressibility such as p-modes (rapid and slow 
magnetic waves, when a magnetic field is applied). Classi- 
cal MHD approximations are used (no charge separation, 
non-relativistic motions, etc.). 



2.1 Equations of motion 

The shell is bathed by an axisymmetric dipolar magnetic 
field: 



Bn — B v ■ R 



)s 8 sin 9 _ 



(1) 



where (e r ,e$,e v ) are spherical coordinate unit vectors. B p 
is the amplitude of the magnetic field at the surface poles of 
the star. 

The equations of motion for an incompressible rotating 
fluid can be written as 



— + v ■ Vv + 2f2 x v = 

at 

-Vp/p - V$ off H (V x B) x B + z/Av, 

Mop 

c?B 

— = V x (v x B) + v m AB 
at 



V B 

V v = 



= 0, 

o, 



(2a) 

(2b) 

(2c) 
(2d) 



where <3> c fr = <E> — il 2 r 2 sin 2 0/2 is the effective gravito- 
rotational potential and fio is the permeability of vacuum. 
Here v is the velocity field, B is the total magnetic field 
(dipole field plus perturbation), and v m = l/(/ioo") where a 
is the fluid's electrical conductivity. 

A non-dimensional form of the previous equations can 
be obtained by introducing the parameters 

V A = Bvl-JWP, Vn = 2QR, Le = V A /Vn, 

E = v/(RVn), E m = is m /(RVn), 

where Va is the Alfven speed, Vn is rotational speed, E is 
the Ekman number and E m the magnetic Ekman number. 



© 2011 RAS, MNRAS 000, 1-9 



3 



Following Jault (2008), we introduce Le, the Lehnert num- 
ber in reference to the work of Lehnert (1954). As shown by 
its definition, this number measures the ratio of the Alfven 
speed to the rotation speed. 

By linearizing magneto-hydrodynamic equations (2) to 
the kinetic u and magnetic b perturbations, one can study 
infinitesimal perturbations from the equilibrium. We assume 
these perturbations have a time dependence of the form e At , 
where A = r + iui (r is the damping rate, w the pulsation 
and i 2 = —1), and we use the following non-dimensional 
variables: 

v -> Vnu(r)e x \ 

B -> BpB(r) + B p b(r)e xt , (4) 

P -> Po, 

where u(r) and b(r) are now first order non-dimensional 
quantities. Therefore, Eqs. (2) reduce to the following set of 
equations: 

AV x u + V x (e z x u) = 

Le 2 V x ((V x b) x B) + EV x Au, (5a) 

Ab = V x (u x B) + E m Ah, (5b) 

V • u = 0, (5c) 

V ■ b = 0, (5d) 

where B denotes the permanent dipolar magnetic field and 
we used V x B = 0. Here we take the curl of momentum 
equation in order to eliminate pressure term. 

2.2 Boundary conditions 

Six boundary conditions are required in order to solve Eqs. 
(5) uniquely. On the inner boundary r = r/R (77 < 1), the 
magnetic field perturbations have only tangent components 
because the core is assumed to be infinitely conducting. On 
the surface r = R, the total magnetic field vector matches 
the external field which is dipolar potential, as there are no 
currents in the external vacuum. As for the velocity field, wc 
may either use stress-free or no-slip boundary conditions. 

As for the magnetic field, different conditions apply for 
the inner and outer boundaries. On the interior, the pertur- 
bation to the electric field is perpendicular to the conducting 
core, and the perturbation to the magnetic field is tangent. 
This gives the following three equations: 



b r 



0, 



(6a) 
(6b) 

(6c) 



— -i-(rbe) = -u B r , 
r dr 

Em 9 . 

r or 

But actually only (6b) or (6c) together with (6a) is needed 
(see Reese et al. 2004). 

The magnetic field perturbations outside the star (r > 
R) are derived from a potential that does not diverge at 
infinity: 

b oxt = V0, (7) 



The boundary conditions at the surface of the star are just 
the continuity of the magnetic field there. These conditions 
are easily expressed after expansion of the fields in spherical 
harmonics (see appendix). 

Eq. (5), together with boundary conditions, defines a 
generalized eigenvalue problem, where A is the eigenvalue 
and (u, b) is the eigenvector which can be computed nu- 
merically. 



2.3 The case of neutron stars 

In neutron stars, typical values for the various physical quan- 
tities are 



B p ~ 10 12 G, p ~ 10 17 kg m" 3 , R ~ 10 km, (8a) 
Q = 2nv s ~ 1900 rad s" 1 (8b) 
a ~ 5.2 x 10 17 ohm" 1 !!!" 1 , ^~0.32m 2 s"\ (8c) 



v m ~ 1.5 x 10 12 m 2 s 1 . 



(8d) 



Here, the values for a and v are given for temperature 
T ~ 10 s K (Baym et al. 1969; Flowers & Itoh 1976, 1979). 
Therefore, Va ~ 2.8 x 10 2 m s~\ Va ~ 3.8 xlO'ms" 1 , so 
that 



Le ~ 7.5 x 10" 



£~8.4x 10" 



E„ 



4 x 10" 



We note that the magnetic Prandtl number v/v m = 
E/E m ~ 10 11 is extremely large. This means that the dif- 
fusion of magnetic perturbations is a negligible source of 
dissipation. As a consequence, we may simplify the set of 
equations (5) by setting E m = 0. In this case the magnetic 
perturbation is readily given by the fluid flow, namely 

\b = V x (u x B) 
Hence, the set of equations reduces to 



AV x u + V x (e z x u) = 

^- V x [(V x V x (u x Bj) x b] + EV x Au, (9a) 

V ■ u = 0, (9b) 

This system is completed by boundary conditions on the ve- 
locity solely. Indeed, the magnetic field is completely frozen 
in the fluid and the magnetic perturbations just repre- 
sent the oscillations of the dipole field lines. Interestingly, 
boundary conditions on the magnetic field (6a) show that if 
E m — then ue = u v = 0. This means that on the inner 
core boundary, no slip boundary conditions can be applied. 
This is expected since the core is assumed at rest and no 
motion of the field lines is allowed. We also assume no-slip 
boundary conditions (u = 0) on the outer boundary to take 
the crust into account. 

The small value of the Lehnert number suggest that 
the coupling between the r-modes and the magnetic field is 
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quite weak. We readily see from the simplified perturbations 
equations (9) that the influence of Lorentz force will be no- 
ticeable compared to the viscous one, if Le > \[E. Noting 
that the Lorentz operator and the viscous operator are both 
of second order, we observe that no length scale comes into 
this inequality. This means that it is valid both in the Ekman 
boundary layers and in the bulk of the layer. From the num- 
bers that characterize neutron stars it is obvious that the 
inequality is met. This means that the magnetic field influ- 
ences the flow more than the viscosity and possibly changes 
the instability of the r-modes. A more detailed calculation 
is therefore necessary to assess the effect. Actually, we shall 
see that the previous inequality is not stringent enough and 
a larger Lehnert number is necessary to affect the instability. 
We now turn to a numerical study of this problem. 



3 NUMERICAL METHOD 

3.1 Spherical harmonic projection 

To solve the eigenvalue problem expressed by Eq. (5), we 
project the set of equations on the spherical harmonics in 
a similar way as in Rieutord (1987, 1991). We expand the 
perturbed velocity and magnetic fields into poloidal and 
toroidal components: 



V = E E uUr)B? +vL(r)ST + wi(r)TT, (10a) 

1=0 m = -l 
oo I 

b = E E a Ur)-RT + bt(r)S7 + cL(r)TT, (10b) 



1=0 m=-( 



where the radial functions {it™,i>m} {{dm,b l m }) and {w l m } 
({ c m}) are the poloidal and toroidal parts of the velocity 
(magnetic) fields, respectively. R™, S™, and T™ are the 
vectorial spherical harmonics 



Rj' 



e. e r , S e = r\/Y e , T, =rVxR, 



(11) 



The harmonic decomposition of Eqs. (5), (6a) and (7) is 
given in appendix. 

Equation (5) reduces to a generalized eigenvalue prob- 
lem 



[A]x = A[B]x, 



(12) 



where [A] and [B] are differential operators with respect to 
the r variable only. The eigenvector associated with A can 
be written as 



Xjm(0 



t+1 (r) 



(13) 



a 



+ E=10"" 
X E = 3*10" 
E = 7*1(T 



□ E = 

nE= 





E=10 1+ i 
E=3*l(Tl>K = 
E=7*10~ 7 = 

E=10~ e o j 
E=3»10" S A = 
E=7*l(T e n z 

E=10" 5 x E 



Figure 1. Top : Absolute value of the damping rate of the m = 2 
r-modc as a function of the Lehnert number for various Ekman 
numbers in [10 -7 , 10 -4 ]. Bottom: The curves TE(Le) are rescaled 
by the critical Lehnert number. The solid line shows the fit \t\ = 
ax 4 +bx 2 +2A5 with a = 1.67, b = -1.63. In both panels j) = 0.35. 



where I is running from m to oo. 



3.2 Classification and symmetries 

Because of the axisymmetry of background fields, the dif- 
ferent m's of the spherical harmonic decomposition are not 
coupled. An additional simplification comes from the sym- 
metry with respect to equator which permits the separation 
of symmetric and antisymmetric modes. A more detailed dis- 
cussion of these points can be found in Reese et al. (2004) 
in the case of pure Alfven modes. In the present case, the 
Coriolis acceleration removes the symmetry m/ — m, which 
exists for pure Alfven modes. Moreover, in the axisymmet- 
ric case, poloidal and toroidal components of the fields are 
coupled. 



3.3 Numerical aspects 

The numerical resolution of the equations is described in 
Reese et al. (2004). We just recall here that the equations 



v m ,w 



. a 



arc 



governing the radial functions u\ 
discretized using the Gauss-Lobatto grid and the resulting 
eigenvalue problem is solved either with a QZ method or 
with the Arnoldi-Chebyshev algorithm according to whether 
we solve for the complete spectrum or a few eigenvalues. 



© 2011 RAS, MNRAS 000, 1-9 



5 



E = 


10" 


+ 


E = 3 


10" 




E=7 


10" 




E- 


10" 




E-li 


10" 




E-'i 


10" 


□ 


E- 


10" 


n 



"A A A A " Aa A a V, 
* • • •••■•! * 



Le/VLe, 




Figure 2. Dependence of the frequency of the m = 2 r-mode as 
a function of the Lehnert number for various Ekman numbers in 
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Figure 3. The critical Lehnert number as a function of the Ek- 
man number for various aspect ratios rj. The straight line shows 
the _E 1 / 4 -law. 



4 i?-MODES-ALFVEN WAVES COUPLING 

r-modes are a subclass of inertial modes which are purely 
toroidal. In the case of a non-magnetic and inviscid incom- 
pressible fluid, exact analytical solutions exist (Rieutord & 
Valdettaro 1997) and the associated velocity field can be 
written as 



u = aQR{r/R) e Ti 



(14) 



where a is an arbitrary constant. The mode's frequency in 
the frame corotating with the fluid is given by 



u = -20. /{m + 1). 



(15) 



Here, we shall focus on the m = 2 r-mode, which is the most 
unstable when coupled to gravitational radiation, and track 
its eigenfrequency and damping rate as the magnetic field is 
increased. 



4.1 The critical Lehnert number 

Since we assumed an infinite magnetic Prandtl number, we 
noted that the boundary conditions on the inner core bound- 
ary are of no-slip type. This means that for Le — > the 
damping rate follows the law derived in Rieutord (2001), 
namely 



35 1 + r/ 
' 2 7 /2 i _ r , 



-I2VE, with I 2 ~ 0.804 



If 77 = 0.35, t = -2.494-/E. 

In Fig. 1, we plot the damping rates for various Ek- 
man and Lehnert numbers. The curves show that when the 
Lehnert number is increased, that is when the magnetic field 
is increased, the (absolute value of the) damping rate first 
decreases before rising rapidly. We define a critical Lehnert 
number, Le c , such that \r(Le > Le c )\ > |i"(Le = 0)|. When 
r is plotted as a function of the rescaled Lehnert number, 
namely Le/Le c , all the curves superpose. 

For later use, we fit this curve with a simple polynomial 
in x = (Le/Le c ) 2/3 , namely \r\ = ax 4 + bx 2 + 2.45 with 
a = 1.67, & = —1.63. The precise shape of the fit is not 
crucial as long as the values of the minimum and the growth 
beyond it are respected. 

In Fig. 2 we show the variations of the frequency of 
the r-mode with both the Ekman number and the Lehn- 
ert number. The curves have been rescaled so as to show a 
minimized dependence to the parameters. This plot suggests 
some scaling laws for w — u>o, notably that lo ~ cjo + oJiLe 4 , 
ui\ > for weak magnetic fields. The asymptotic analysis 
of the behaviour of uj is quite cumbersome and beyond the 
scope of this paper. 

We show in Fig. 3 that the critical Lehnert number 
varies quite closely as E 1 ^. Actually, a good fit is Le c = 
1.78-E 1 / 4 . Magnetic fields will therefore modify the dynamics 
of the oscillations when Le ^ Le c . 

The E 1 / 4 scaling law of the critical Lehnert number be- 
yond which magnetic field is influential is rather surprising 
in view of the discussion of sect. 2.3. We may understand 
this scaling if we consider the nature of r-modes when per- 
turbed by a dipolar magnetic field. The magnetic and veloc- 
ity perturbations, as long as their gradients, are displayed in 
Fig. 4. It is clear that the magnetic perturbation is concen- 
trated along the shear layer emitted at the critical latitude, 
namely arcsin(l/3) = 19.5°. Such shear layers are a feature 
of viscous inertial modes which is due to a singularity of the 
boundary layer at this latitude (e.g. Rieutord & Valdettaro 
2010). 

By analyzing the scale of this layer, as shown in ap- 
pendix, we find that when Le S> VE, its thickness scales 
like Le/E 1 / 4 , showing that indeed when Le reaches E 1 ^ 4 , 
the interaction between the magnetic field and the r-mode 
becomes strong, presumably emphasizing a resonance be- 
tween the oscillating shear layer and an Alfven wave. From 
the diminishing of the damping rate, we conclude that the 
dissipative layers slightly thicken, which is understandable 
since the magnetic field tends to oppose to shear. 

When the Lehnert number is much larger than its crit- 
ical value the r-mode is completely destroyed, leaving the 
place to small scale (very dissipative) Alfven waves (see 
Fig. 5). 
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Figure 4. Top : In a meridional plane of the shell, the kinetic 
energy (right) and viscous dissipation (left) of the m=2 r-mode 
when a weak magnetic field is present (small Lehnert number). 
Bottom: a view of the associated magnetic perturbation. On right 
the magnetic energy associated with the oscillation. Left: the as- 
sociated current. Here E=10 — 6 , Le=10 — 3 and r\ = 0.35. 




Figure 5. As in figure 4 but with a critical Lehnert number. Note 
how the dipolar magnetic field now features the various part of 
the oscillation. 



4.2 Influence of the magnetic field on the 
instability 

Let us now turn to the physical consequences of the inequal- 
ity 

Le ^ 1.78£ 1/4 . 

It may be transformed into a constraint on the magnetic 
field, namely 

B > Bcrit = 1.78^Jtip2QRE 1/i . 

We now introduce the scaled angular velocity f2„ = fl/ 'y/nGp 
and the temperature dependence of the kinematic viscos- 
ity as in Bildsten & Ushomirsky (2000), namely v = 
1.8//T| m 2 s _1 , where / is a parameter of order unity and 
T n = T/10". We thus find 

O 3 / 4 

S crit =7.110 14 ^ 7F G. 

This critical magnetic field is very high, showing that only 
hot and slowly rotating neutron stars might be affected by 
the magnetic field. To visualize this effect, it is useful to 
consider the window of instability in the plane (T, fi„) as in 
Rieutord (2001) or in Bildsten & Ushomirsky (2000). 

In order to calculate the boundaries of the window, we 
need to approximate the curves in Fig. 1 with some analyt- 
ical function. We thus use the polynomial derived preced- 
ingly. 

In the plane (T, Q*), the window boundary gives the 
critical angular velocity beyond which the r-mode instability 
exists. We therefore need to solve 

Jgw + jBulk + 7mi = 

where we take the growth rate of the mode due to grav- 
itational radiation j gw — 0.306 s^ 1 Q,^ from (Lindblom 
et al. 1998) and the damping by bulk viscosity jBuik = 
-2.2 10 -12 s _1 T 9 6 fi 2 from Lindblom et al. (1999). These 
expressions of ~/ gw and "/Bulk have been derived using an 
n=l-polytropic model for the neutron star, however we ex- 
pect that the density distribution influences the damping 
rate of the mode with a factor of order unity, thus not chang- 
ing the orders of magnitude of the magnetic fields. We can 
thus estimate in a simple manner the critical rotation rate 
for various values of the temperature and magnetic field. 

As expected from the value of the critical magnetic field, 
we see that there is a significant reduction of the instability 
window only when the field exceeds 10 14 G, which is rather 
high. However, we note that for lower values (e.g. 10 13 to 
4 10 13 G), the magnetic field slightly widens the instability 
window. This is the consequence of the slight reduction of 
the damping rate around the critical Lehnert number that 
appears in Fig. 1. 
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Figure 6. The critical angular velocity, normalized by Qk = 
%^\J-kG~p for various values of the magnetic field. The solid line 
shows the case with no magnetic field. 



& Drago (2010) in the context of accreting neutron stars 
like low-mass X ray binaries. Neglecting viscosity, these au- 
thors show that magnetic fields above 10 12 G can reduce 
the r-mode instability. From these results and ours, we es- 
timate that accounting of viscosity will require fields ten 
times stronger, namely above 10 G. However, the details 
of the interaction between r-modes and a toroidal field in 
a spherical fluid layer are poorly known: only a few works 
(e.g. Schmitt 2010, and references therein) have investigated 
this question in the context of the magnetohydrodynamics 
of the liquid core of the Earth, where the magnetic Prandtl 
number is very low. The case relevant to neutron stars still 
deserves detailed investigations. 



The numerical calculations have been carried out on 
the NEC SX8 of the 'Institut du Developpement et des 
Ressources en Informatique Scientifique' (IDRIS) and on the 
CalMip machine of the 'Centre Interuniversitaire de Cal- 
cul de Toulouse' (CICT) which are both gratefully acknowl- 
edged. 



5 CONCLUSIONS 



In this paper we investigated the dynamics of the m = 2 
r-mode in a spherical shell when the fluid is bathed by a 
dipolar magnetic field. The main theoretical result of this 
work is that magnetic fields have a negligible influence on 
the r-mode unless they exceed a critical value which is such 
that the Lehnert number is of the order of the 1/4-power of 
the Ekman number. A physical interpretation of this scaling 
law may be obtained by observing that at this critical value, 
the wavelength of Alfven waves is similar to RE 1 / 4 which 
is precisely one of the typical thicknesses of shear layers of 
inertial modes (e.g. Rieutord et al. 2001). 

We also have shown that when the magnetic field is 
slightly below the critical value, the damping rate of the 
mode is slightly reduced because of the thickening of the 
various layers due to the frozen field limit. 

Back to neutron stars and the famous instability of in- 
ertial r-modes when coupled to gravitational radiation, we 
conclude that, as far as dipolar fields are concerned, only 
fields over 10 14 G may seriously affect the instability. Thus, 
with a quite different approach, which includes the specific 
shear layers of r-modes, we can confirm the result of Lee 
(2005), obtained with singular perturbations, that for typi- 
cal values of 10 12 G, magnetic fields are unimportant. 

This result does not therefore invalidate the scenario 
of Rezzolla et al. (2001b), which suggests that the nonlinear 
development of the r-mode instability may generate a strong 
toroidal field from a pre-existing poloidal field. This mech- 
anism was recently investigated in some details by Cuofano 



APPENDIX A: SPHERICAL HARMONIC 
EXPANSION OF THE MHD EQUATIONS 

In this section , we expand each term of Eq. (5) in the spher- 
ical harmonics in base of (R™ , S™ ,T™). The different parts 
of the equation have been published elsewhere: The Lorentz 
force and the induction equation may be found in Rincon 
& Rieutord (2003) while the projection of Coriolis accelera- 
tion or its curl are in Rieutord (1987); Rieutord & Valdettaro 
(1997). For completeness we give here the result of casting 
the projection of these forces into a single equation. The 
equation of vorticity yields: 



Aw 



im i i ;_i d ( u m 
__ w +Al _ lX 



A -1-2 O i 

-A M x — \x 



dx 



1+3 1+1 



+B e e+1 d r (rc e + 1 ) - 



-rA e (ra m ) 



Be 
r 



[(£ + Vatic*' 1 - ea e e+1 c e +^ | + EA e wi 



(Al) 
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XA e (ru e m ) - A e (ru e m ) + 



1(1 + 1) 



= Le 2 £(l + 1) |^ [iAtiAe-^ra*- 1 ) 

-{I + l)^ +1 A m (ra^ 1 ) + imct/r] 
-B e imc m /r 

+B r ^A(_ 1 T>(-iaf n 1 + A e e+1 V e+1 a e „ 



— I r w, 
dr 



+ 



im Id ret 



1(1 + 1) r dr 2 

+B' r [A t t _ 1 A e -i{ra e - 1 ) + A t t+1 A i+ i(ra t + 1 ) + 

im 1 drcil 1 < 
+ SA^A^rw™) 



1(1 + 1) r dr 
While the induction equation gives: 



(A2) 



A-i = (i - l)a{_!, A'+i = + 2 W+i 
7,* = ((*(* + 1) + 1(1 + 1) - 2)/2fc(fc + l))a? 

= A'- 1 = (/ 2 - l)a{_! 

All the differential equations need to be completed by 
boundary conditions; for the velocity field, no-slip conditions 
read: 



U m = —r = W m = 

ar 

on both boundaries. These conditions are sufficient when 
E m = 0. 



APPENDIX B: BOUNDARY LAYER ANALYSIS 
FOR THE CRITICAL LEHNERT NUMBER 



\ra„ 



Ac„ 



imB r w e m + Be \ (l + X)ci S i_ x u m 1 — la£ +1 u^ 
B r [bI^v 1 ^ 1 + B l l+1 v l +^ + E m A t (ra m ) 



r r 



wt + BA^ 
dr 

i 



^ T) (e(e + i)B g + (rB r y) v ^ r 

im dv l m (rB e )' im e 

-£>T ; 1 „x ^n 



£(£+!) dr r £(£ + 1 ) 

R dui m , p A j 

and the mass and flux conservation read: 



1 d -(r 2 u e m ), (V-v = 0), (A3) 



+ l)r dr 

d (r 2 ai), (V-b = 0), (A4) 



m <?(<?+ l)rdr 
In these equations we introduced the following operators 



Let us consider cartesian coordinates adapted to the geom- 
etry of the shear layer as shown in Fig. 4 bottom (as ma- 
terialized by the region with high b 2 ). z is the coordinate 
normal to the shear layer assumed to contain the rapid spa- 
tial variations. 

The dynamics inside this layer verifies 



Xu + uje z A u = — Vp + 



he 



2 ' du^ 

V x — ] x e z + EAu 

oz 



where we assumed that the background magnetic field is 
along e z and that the shear layer makes an angle arcsinw 
with the rotation axis (hence the factor u in front of Coriolis 
acceleration). After taking the curl of this equation, consid- 
ering only the rapid variations along e z , and taking its x 
and y components, we find 

A(A + iu))d z u = (he 2 - \E)dlu 

with u = u x + iu y . 

Noting that for r-mode A = — iu> + r \/E, we find that 
if u oc exp(— z/zo) then 



f d 2 1(1+1) 



r dr z 



-r, f d , f d 3 
V t = --j-rAe = --Ter- 
ror r dr 6 

and coupling coefficients: 



'+ 1) d f 



r dr r 



I 2 - m' 2 



(2/-l)(2Z + l) 



^~ ( 1 + l )^i7I 

This result shows that there is a change of scale in the flow 
when the Lehnert number reaches values 0(E 1 ^ 4 ). 
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